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Abstract. We study the long time behavior (homogenization) of a diffusion in random medium with time and space 
dependent coefficients. The diffusion coefficient may degenerate. In Stochastic Process. Appl. (2007) (to appear), an 
invariance principle is proved for the critical rescaling of the diffusion. Here, we generalize this approach to diffusions 
whose space-time scaling differs from the critical one. 

Resume. Nous etudions le comportement asymptotique (homogeneisation) d'une diffusion en milieu aleatoire avec 
des coefficients dependant du temps et de I'espace, pour laquelle le coefficient de diffusion pent degenerer. Dans 
Stochastic Process. Appl. (2007) (to appear), un principe d'invariance est etabli pour le changement d'echelle critique 
de la diffusion. Ici, une generalisation de cette approche est proposee pour differents changements d'echelle possibles. 



1. Introduction 

Wc aim at studying the long time behavior, as e — > 0, of the diffusion process with random coefficients (the 
parameter lu below stands for this randomness) defined by 



and then at identifying the limit of the solutions of the random parabolic differential equations (PDE) 



dtZe,ij{x,t) = ifrace 



+ [e-/^c-\-d](^-^,^,u?jz,^^{x,t), {x,t)eR'' X 



(2) 



with initial condition Zi,^^(x,0) = f{x). a, (3 are two strictly positive parameters. The coefficients a,b,cr,c,d 
are stationary ergodic random fields with respect to space and time variables. We shall see that there are 
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different possible limits depending on a — 2/3, a < 2/3 or a > 2/3. More precisely, we suppose that a = aa* 
and the generator of the diffusion process could be rewritten in divergence form as 

L-^-=Q)e2^(-/^'^-)div,(^e-2^(-/^''-)[a + i7](^i^,^,c.)v.). (3) 

Here V and H are also stationary random fields and H is antisymmetric. Assumptions will be stated 
rigorously in the next section. We will prove that, in probability with respect to w, 

lim Zs,ujix,t) =z{t,x), (4) 

e— »0 

where z is the solution of a deterministic equation of the type 

dtz{x, t) = iv&ce[AAj,r^z] (x, t) + C • ^^zix, t) + Uz{x, t) (5) 

with initial condition 'z{x, 0) = f{x). A, C, U arc deterministic coefficients, the so-called effective coefficients, 
and depend on the considered case. As in the previous paper [18], we do not assume any non-degeneracy of 
the diffusion matrix a (several examples are introduced in [18]). 

Roughly speaking, the case a = 2/3 corresponds to the critical scaling for the diffusion and space and 
time variables have to be homogenized simultaneously. This job was carried out in [18]. In case a < 2/3, the 
space variables are moving faster than the time variable so that the homogenization procedure has to be 
performed first in space and then in time, and vice versa in the case a > 2/3. 

Both this paper and [18] follow a series of works on this topic. Let us sum up the methodological ap- 
proach of this issue. Briefiy, the time dependence of the process brings about a strong time degeneracy of 
the underlying Dirichlet form, which satisfies no sector condition, even weak. To face such a degeneracy, 
analytical methods, more precisely compactness methods, carried through the homogenization procedure 
for a uniformly elliptic diffusion matrix a in periodic media [1, 20] or in random media [4]. In [8], critical 
scaling is considered {a — 2 and P — 1) and an annealed invariance principle is proved for the diffusion (1) 
by means of probabilistic tools under the assumptions V = and a = Id. This result is extended in [10] to 
the case when the stream matrix H satisfies a certain integrability condition of the spatial energy spectrum 
instead of boundedness of the coefficients. A quenched version of the invariance principle is stated in [5] for 
an arbitrary time and space dependent diffusion coefficient provided that it satisfies a strong uniform non- 
degeneracy assumption. The particularity of these works lies in their intensive use of regularity properties 
of the heat kernel to deal with the time degeneracy of the Dirichlet form. 

In this paper, we additionally consider possibly degenerate diffusion coefficients a. This prevents us from 
using both compactness methods (the Poincare inequality is lacking) and regularity properties of the heat 
kernel (no uniform ellipticity or even hypo-ellipticity of the matrix a is assumed). The main tools of the 
proof are the control of the matrix a by a time independent one 5 and a commutativity argument of the 
unbounded operators associated to d and to the time evolution. 

This article is structured as follows. In Section 2, we introduce notations and present our results. In 
Section 3, we set out the main properties of the involved stochastic processes. Section 4 is devoted to 
constructing the correctors. We explain how to get round the lack of regularity of the correctors in Section 
5 to establish the Ito formula. Then in Section 6, an ergodic theorem is proved, which allows us to carry 
through the homogenization procedure (Section 7). The tightness of X'^'^ is studied in Section 8 by means 
of the Garsia-Rodemich-Rumscy inequality. The limit PDE 5 is identified in Section 9 with the help of the 
Girsanov transform. 

2. Notations, setup and main results 

The setup remains the same as in [18]. It is reminded for the reader's convenience. 

Definition 2.1. Let {n,Q,fi) be a probability space and {rt^x', {t,x) G M x M**} a stochastically continuous 
group of measure preserving transformations acting ergodically on f2: 
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(1) Vyleg,V(t,x)eMxR^ ^(Tt,,A) = 

(2) If for any {t,x) € R x E'', Tt^^A^A then //(A) = or I, 

(3) For any measurable function g on {n,Q,iJ.), the function {t,x,uj) i— >g(Tt^a;Ci;) is measurable on (R X M'^ X 
i7,S(Rx R'^)®^;). 

In what follows wc will use the bold type to denote a function g from H into R (or more generally into R", 
n > 1) and the unbold type g{t,x,uj) to denote the associated representation mapping (t,x,u!) i-^ g{Tt^xi^). 
The space of square integrable functions on {n,0,ij.) is denoted by L^(i7), the usual norm by | • I2 and the 
corresponding inner product by (•,-)2- Then, the operators on L'^{f2) defined by 2:g(a;) = g(Tt^j.w) form 
a strongly continuous group of unitary maps in L'^{Q). Each function g in L'^{f2) defines in this way a 
stationary crgodic random field on R'^^^. The group possesses d + 1 generators defined for i = 1, . . . , c?, by 

d d 
Dif ^ -^To^xi\{t,x)=o and Dti = —Ttfii\(^t,x}=o, 

which are closed and densely defined. Denote by C the dense subset of L^{fl) defined by 

C = Span{f *v3;f eL^^^^^^gC-oo^jgd+i^l^ withf*v5(w)= / f(Tt,,,tj)^(t,x) dtdx, 

where C^iW^^"^) is the set of smooth functions on R'^+^ with a compact support. Remark that C C Dom(_Dj) 
and -D,;(f * 1^9) = — f * This last quantity is also equal to ZJ^f * if f e Dom(_Di). 

Consider now measurable functions cr, or, H : i? ^ R''^'' and V,c,d:J?^R. Assume that H is anti- 
symmetric. Define a = crcr* and a = aa* . The function V does not depend on time, that means \/t G R, 

rt,oV = v. 

Assumption 2.2 (Regularity of the coefficients) . 



• Assume that\/i,j,k, I = 1, ... ,d, aij,a.ij,Y,llij,Diaij, and Diaij GDom(Dk). 



"^^Ei., A(e-2^ff,jfj) for some function { e L°°{f2;R'^) 



• Define, for i = 1, . . . ,d, 



and assume that the applications {t,x)^hi{t^x,U}), {t^x)^~* c{t,x,uj), {t,x)^ a{t,x,uj), x^-^ DV{x,u}) are 
globally Lipschitz (for each fixed co € f2) and the application (t, x) d{t,x,uj) is continuous. Moreover, 
the coefficients cr, a, b, or^ V, H, c, d are bounded by a constant K. (In particular, this ensures existence 
and uniqueness of a global solution of SDE (1).) 

Here is the main assumption of this paper: 

Assumption 2.3 (Control of the coefficients). 

• or does not depend on time (i.e. \/t G M, T^or = orj. As a consequence, the matrix a does not depend on 
time either. 



H, a, c G Dom(Df) and there exist five positive constants to(> 0), A/, Cf^ , C2 , such that, fi a.s., 

msL<a< Ma, (7) 
|H| < Cf S, I AH| < Cf S and | Aa| < C^a, (8) 
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where |A| stands for the symmetric positive square root of A, i.e. |A| =\/ AA* . 

For instance, if the matrix a is uniformly elliptic and bounded, or can be chosen as equal to the identity 
matrix Id and then (8) ^ H, AH and Aa G 

Let us now set out the crgodic properties of this framework: 

Assumption 2-4 (Ergodicity). Let us consider the operator S = (l/2)e^^^^ - j Z)i(e~^^aijA) 
domain C. From Assumption 2.2, we can consider its Friedrich extension (see [6j, Chapter 3, Section 3), 
which is still denoted S. Assume that each function f G Dom(S) satisfying Sf = must be fi almost surely 
equal to some function that is invariant under space translations. 

The reader is referred to [18] for examples in which Assumptions 2.2, 2.3 and 2.4 are satisfied. 

Even if it means adding to V a constant (and this does not change the drift b, see (6)), we assume that 
j e~^^ d/i = 1. Thus we can define a new probability measure on fl by 

d7r(cj)=e-2V(-)d^(^). 

We now consider a standard d-dimensional Brownian motion defined on a probability space P) 
and the diffusion in random medium given by (1). Under Assumptions 2.2, 2.3 and 2.4, we claim: 

Theorem 2.5. The law of the process X*^'" converges weakly in ^-probability in C([0,T];R'') to the law of 
a Brownian motion with a certain covariance matrix (see (50)). 

This result puts us in position to describe the long time behavior of PDE (2): 

Theorem 2.6. For any bounded continuous function /:]R'^'^M, the solution z^i^{x,t) of PDE (2) with 
initial condition Z£.i^(a;,0) = f{x) converges in ^-probability towards z{x,t), where z is the unique viscosity 
solution of the deterministic PDE (5) with initial condition z(.t,0) = f{x) (see (49) and (50) for a description 
of A, C, U ). More precisely, for any (x, t) G M'' x and S > 

\zs,^{x,t) - z{x,t)\ >5})^0 

as s tends to 0. Furthermore, the coefficients A,C,U only depend on the case a — 2/3 < 0, a — 2/3 > and 
a -2/3 = 0. 

The reader can refer to [2] for a description of the viscosity solution theory of PDEs. 

Remark 2. 7. Let us be more explicit about the last statement of Theorem 2. 6. Because of the possible 
degeneracies of the diffusion matrix a, the coefficients A, C, U are defined by limits involving solutions of 
parameterized PDEs (see Section 9). However, the uniformly elliptic setting provides us with exact problems 
(or PDEs) to compute A, C,U . The reader is referred to [1], Chapter 2, Section 1, or [16, 20] for the periodic 
case and to [4] for the random case. 

To get an idea of the scaling effect, let us focus on the 1-dimensional periodic case. In this setting, the 
matrix A is given by (see [20], Proposition 2.2) 

A= / a{t,x){Dv{t,x)-\-l)dxdt, (9) 

Jt£lO,l] Jx£lO,l] 

where the function v solves one of the following problems, depending on the values of the parameters a and 
/3. 

• Case a — 2/3 < 0.' V is the unique solution of the parameter dependent elliptic problem: 
^D,{a{t,x){D,v{t,x) + l))^0, 

v(t, •) is 1 -periodic, / v(t,x)dx^O, tG[0, 1]. 

Jx£[0,l] 
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It is readily seen that Dx'v{t,x) ~ C{t)/a.{t,x) — 1 where C is a function of t. Since v is x-periodic, 
/xe[o 1] -Da;V(<, x) dx = 0. Consequently, it comes out that C{t) ~ (/^.^[o i] g,(t x) ^^)~^ that, from (9), we 
obtain explicitly the effective diffusivity 

f ( f — ^ dx\ dt. 

Jte[o,i] \Jx€lo,i] a.[t,x) ) 

• Case a — 2/3 > 0; V is the unique solution of the elliptic problem: 
-D,{k{x){DM^) + l)) = 0, 



V is 1-periodic, / v(a;) dx ~ 0, 
'xe[o,i] 

where a(x) = J^^^^ ^a(t,x) dt. As in the previous case, solving this problem raises no particular difficulty 
and we get Dxv{x) =C/a(x) — 1 where C is a constant, which exactly matches (/^^jq g^ky*^^)"^ thanks 
to the periodicity of w again. Then formula (9) reads 

A~{ I - dx ' 

'xe[o,i] /t6[o,i] 

Note that, in comparison with the previous case, the role of time and space variables have been, in a way, 
inverted. 

• Case a — 2/3 = 0; V is the unique solution of the parabolic problem: 

DMt, x) - D^mt, x){D^v{t, x) + 1)) = 0, 

V is periodic in time and space, ft j,^[q v(t, x) di dx = 0. 

Unfortunately, computations in the critical case a = 2(3 are trickier and do not enable us to give an explicit 
formula for the effective diffusivity. 

From now on, in order to avoid heavy notations, the superscript lo of X"^'^ is omitted when there is no 
possible confusion. So we write X^ instead of X^''^. Moreover, except where otherwise stated, the process 
X^ starts at time s = {) (sec (1)). 



3. Environment as seen from the particle 

Let us denote by X^ the process driven by the following Ito equation 

Xl = [ b{e^f^-°'r,xl,Lo)dr+ f a(e^^-°'r,xl,uj)dBr. (10) 
Jo Jo 

From law uniqueness for solutions of Ito equations with Lipschitz coefficients, the processes 
{e^Xf/^2i3;t > 0} and {X^;t > 0}, which both start from 0, have the same law. 

We now focus on the environment as seen from the particle: this is a process on the probability space f2 
defined for each e > by 

where the process X"^ starts from S M'^. This section is devoted to proving that it defines a valued Markov 
process, which admits n as (not necessarily unique) invariant measure. As a consequence, the associated 
semigroup (see [3] for a thorough study on semigroups) can be extended to a contraction semigroup on 
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U'{Q). for l<p< oo. Let us additionally mention that, with the help of the Ito formula, the generator of 
this process is easily identified on C. It coincides with the operator L + e^'^^'^L't, where L is defined on C by 

2V d 

L = — 5] A(e-2^[a + H],,i?,). (12) 

ij = l 

Note that C need not be a core for this operator because of the possible degeneracies of the matrix a. 

All these statements are well known in the case of a uniform elliptic diffusion coefficient a for time 
independent (see [9, 13, 14, 15]) or time dependent (see [8, 10] and references therein) coefficients. In what 
follows, we will see that the degenerate case boils down to the uniform elliptic case by means of vanishing 
viscosity methods. 

Let us consider a ((i+ l)-dimensional Brownian motion B' independent of B. Up to the end of this section, 
the couple (a, X), a G M, X g M'' stands for a. (d+ l)-dimensional vector and the variables m, v denote (d+ 1)- 
dimensional vectors. Let us define the [d + l)-dimensional process X*^'" as the solution of the following Ito 
equation 

dX,"'" = (e'/^-", b - n-^DV{X^'",u)) dt + (0, ct(X,"'", w) dBt) + n~^l^ dB[. (13) 

From Lemma 3.1 below, the R'^+i -valued process X*^'" admits transition densities Pe_„(i, u,?;) with respect 
to the Lebesgue measure on M'^+i. The transition densities satisfy for any t>0 and u,v G K'^+i 

■ exp \ < p,,„ (t, u, v) < ™^ exp ~ J (14) 



for some constant C > that only depends on K^d,n. In particular, we can now easily adapt the proof in 
[9], Section 2.3, and prove that the f2-valued process y/'"(a;) = tx^^'"<^, X^'" starting from at time t = 0, 
is a Markov process whose generator coincides on C with 

e2V JL 

= V ^ A(c-2^[a + H + 72-ild].y.Z?,) + (2n)-iZ?2^e2^-"A, (15) 

and admits tt as invariant measure. 

It now remains to let the parameter n go to oo. For this purpose, define the subspace G{Q) of L°"{Q) as 
follows: f G L°°{il) belongs to C{S1) if and only if it is bounded and for ^ almost every oj G f2, the function 
{t,x) €RxR'^^ f(Tt,:jU;) is continuous on R x W^. Let C{f2) denote the closure of C(i7) in L°°{f2) with 
respect to the usual L°°(I?)-norm. For each function f G C{f2), let us then define 

Pt(f)(c.) =Eo[f(r .,_„,^^.c.)] =E[f(r/(c.))]. 

Obviously, Pt is continuous with respect to the L°°(i7)-norm. Lemma 3.2 below ensures that Pt(C(J7)) C 
C(J7). Then, from the Markov property of the R'^+^-valued process {e^^~°'t,X^) and (18), it is readily seen 
that Ps{Ptf) = Ps+t{f) so that the family {Pt)teM+ defines a semigroup on C(i7). 

Let us now consider a function f G C(]7) and fix a; G i?. From the Lipschitz property of the coefficients 
and the continuity of {t, x) gRxW^ i-^ i{Tt^x(^), classical arguments of SDE theory imply that Eo[f (TjfE."tj)] 
converges to Eo[f (T^2fi-at x^*^)] Soss to oo. Together with the boundedness of f, this ensures that 

tt{Eo[{{tx^-"'Uj)]) converges to 7r(Eo[f (r^s^-at x^^^)]) as n tends to oo. As 7r(Eo[f (Tx^^.^tj)]) =7r(f), it comes 
out 7r(Eo[f (t^23-<:.( y^^)]) = '"'(f)- Hence tt is invariant for the semigroup {Pt)teR+ , which now clearly extends 
to a contraction semigroup on U'{Q,'n) for 1 <p < oo. 

Lemma 3.1. The transition densities with respect to the Lebesgue measure of the process X^'^ satisfy (14)- 
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Proof. Let us define a new {d+ l)-dimcnsional process 

dX^''' = {0,b~n-^DV{X^-",uj))dt+{0,a{X^^'\Lj)dBt)+n-^/^dB',, (16) 
We point out that its generator on C^(R'^+^) can be rewritten in divergence form as 



i.3=0 



(with the convention Oio = aoi = HiQ = H^i ~ 0). From [12], this process admits transition densities with 
respect to the Lebesgue measure satisfying 

■ exp < p,,„(t, X, V) < ™^ exp - (17) 



for some constant C > that only depends on n, K, d. 

Let us now suppose that aU the involved coefficients in (13) or (16) are smooth. For a C^{M.'''~^^) 
function / (of class with bounded derivatives up to order 2), it is well-known that the mappings 
(t,x) G R X R'^+i ^ E^[/(Xt^'")] and {t,x) e M x M''+i ^ E^[/(Xj^'")] ai-e at least of class C^'^iR x K'^+i) 
(see [7]) and are respectively classical solutions of the PDEs dtu = U-'^u and dtu = L^'"'u + e^^~°'doU 
with initial conditions u{0, •) = /. Then, it is readily seen that the difference between these functions 
{t, x) 1-^ E^[f{X^''^)] - E^+e2f!-»^o [/(^f'")] is a classical solution of the PDE dtu = L'-"u + e'^^'^'d^u with 
initial condition u(0, •) =0. Thus, from the comparison theorem (see [17], Theorems 2.4 and 3.1, for a 
probabilistic proof of this fact) the functions coincide on ]R''+^. With density arguments, we establish 
Ex[f{Xl''^)] = E^+j23~Q3,^ [/(X(^'")] for each continuous bounded function /. The lemma then follows from 
(17). 

If the coefficients in (13) or (16) are not smooth, the situation easily comes down to the previous case 
with the help of a regularization procedure. Details are however left to the reader. □ 



Lemma 3.2. For each fixed uj & Q and {s,x) G M x M'^, the process t e M+ ^ (e "(t + s),e ^x + 
Xj^'^='"°'"^"""), starting from x^'^^'''"-"''"''" = 0, and the process t e M+ (£-"(t + s) , Xl;^^) , starting 
from = X, have the same law. As a consequence, for each fixed t € M+ and f G C{Q), the function Pt{f) 
belongs to the space G{f2), and 

Pt (f ) (t,2,-„,^,c.) = [/(£2/5-"(t + s)X+s , ^)] ■ (18) 



Proof. For the sake of clarity and without loss of generality, let us drop the parameter s by choosing £ = 1. 
Fix (s, x) e M X R''. It suffices to check that the process {s + t,x + Xl"'""^), xl"'""^ starting from e R'' at 
time t = 0, and the process {t + s,X^j^^), starting from x S R'^ at time t = Q (that is X^ = x), satisfy 

the same Ito equation, and therefore have the same law by virtue of law uniqueness for solutions of Ito's 
equations. As a consequence, for f S C(/?), the function 

(s, x) ^ E[f (rt(T,,,co))] = E,[/(t + s, X^+,,uj)] 

is continuous from the continuity of f and the Lipschitz properties of the coefficients b and cr. □ 



4. Resolvent equation 



Let us now investigate, for each A > 0, the resolvent equation (ua is the unknown) 



AuA-(L + 6'A)uA = h 



(19) 



680 



R. Rhodes 



for a function h e L?'{Q) and a function 6 = 6{X) of the parameter A. Note that the operator L + 6Dt has 
not been rigorously defined yet but a complete description of all the involved operators and the meaning of 
"solution of (19)" are given thereafter. Actually, because of both time and space degeneracies, the generator 
of the process does not possess enough regularity to work on in a quite general way. However, for a suitable 
right-hand side h, the function \i\ inherits some regularity properties that allow us to carry through the 
study of (19). This argument will be the guiding line of this section. 

The following study is carried out for a quite general strictly positive function 6 of the parameter A. 
However, this result will only be used in Section 7 with 9 = 9{\) = . Roughly speaking, this function 

measures the difference of speed rates between the time and space variations. 



4-.1. Setup 

In [18], it is proved that the imboimded operators S and Dt on L^(I?,7r) (see the definition in Assumption 
2.4) have a common spectral representation. This is due to the time independence of the coefficients a and 
b. More precisely, we can find a spectral resolution of the identity E on the Borelian subsets of M_|_ x M such 
that 



: / xE{Ax,Ay) and —Dt ~ I iyE{dx,dy). 



For any ifi,ip E L^(i7), denote by the measure i?,^,^ = {Eip,xp)2. Let S be the Friedrich extension of 

the operator defined on C by (l/2)e2V^^ . D,{e~'^^ a^jDj). For any ip,xp EC, define 



jR+xR 



and llv'lli = {ip,ip)\^^, \\ip\\i^s = (ViV)!^? the corresponding seminorms, whose kernels both match the 
L^(f2)-sub-spacc of functions that are invariant under space translations (sec Assumption 2.4). By virtue of 
assumption (7), these seminorms are equivalent 

"^ll¥'ll?<il'^llU<^^llvlli- (20) 

Let F (respectively H) be the Hilbert space that is the closure of C in i^(i7) with respect to the inner 
product e (resp. k) defined on C by 

e(¥','0) = (¥',■0)2 + + (Ay, AV')2 

(resp. k(v,'0) = (¥',V')2 + {'P,tp)i)- 

Define the space D as the closure in (L^(i7), | • I2) of the subspace V = {(— S')^/^^?; (p EC}. For any if EC, 
define ^{{-Sy/^(p) = a*D,^if E {L'^{Q)Y. Thanks to the description of the kernel of the semi- norm • h, 
note that <!> is well defined on V. Furthermore |^((— S')^/'^¥')|2 = ^2(c^, 8(^)2. From (20), <P can be extended 
to the whole space B and this extension is a linear isomorphism from D into a closed subset of (i^(i7))''. 
Hence, for each function u e H, we define V^u = <?((— S')^/^u) and this stands, in a way, for the gradient of 
u along the direction a. 

For each f E L'^{ fi) satisfying /u^^r ^-^f.f ('^^j '^y) define 

l|f||-i=/ -E,s{dx,dy). (21) 

JR+ xR 

We point out that ||f < 00 if and only if (see [14], for instance, for further details) there exists C £ M such 
that for any if eC, (f, ^5)2 < C||(p||i. For such a hmction f , ||f also matches the smallest C satisfying this 
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inequality. Remark that ||f || -i < oo implies 7r(f) = 0. Denote by H_i the closure of L^{fi) in H* (topological 
dual of H) with respect to the norm || • 

Let us now focus on the antisymmetric part H. We have 

|(w,Hw)| < (u,|H|u)^/^(w,|H|w)^''^ <Cf (22) 
The second inequality follows from (8) and the first one is a general fact of linear algebra. We deduce 

Vv,t/>GC, {^{liD,^,D,il,)2<C^\mi\Mi- (23) 
For any (p,tp eC, let us define 

{llD,.cp, D,xP)2 = TnihSr^'v, i^Sf'^^l^)- (24) 

Note that, if ip,ip' and {-Sy/'^ip = {-Sf/'^ip' , then - tp') = and, from Assumption 2.4, D,j,ip = 
D^if', in such a way that T//((-S')i/2(p^ (_5;)i/2^) ^ Th((-S')1/V, hS)^^'^tp)- Thus T^- is a well-defined 
antisymmetric bilinear operator on V xV. From (23), it extends to an antisymmetric continuous bilinear 
form Th on D X D. Likewise, with the help of Assumption 2.3, we define the continuous bilinear forms Ta, 
dtTa, dtTn, ^T^, on D x B C L'^ {[2 , ir) x L'^{f2,Tr) as follows: \/(f,ipeC, 



{aD,^, D,i,)2 = T,((-5)i/V, (-<5)i/V), 
{DtaD.,cp, D,'iP)2 - atTa((-^)i/V, {^Sf'^'^l^). 



^ (AHZ?,<^, D,,l,)2 - 9tTH((-5)i/V, {-Sf'^-iP), 



{A^aD^^, D,xP)2 = ^T,((-5)i/V, {-S)'^^ip), 
{AsliD:,V, D^,p)2 = ^Th((-5)1/V, (-^)'/'-0), 



where, for any s G M*, Ag denotes the L^-continuous difference operator (remind of the definition of Tgfi in 
Section 2): 



^(f) = I^MLJ£. (25) 

From Assumption 2.3, the norms of the forms AsTa and A^Th are uniformly bounded with respect to 
s G R* and the forms are weakly convergent respectively towards dtTa and dtTn as s ^ 0. 
Now, denote by H the subspace of EI_i whose elements satisfy the condition: 

3C>0,Vs>0 and Vy-GC, (h,^(p)_i,i < C||v>||i. (26) 

For any h G 7i, the smallest C that satisfies such a condition is denoted ||hj|T. Then H is closed for the 
norm || • \\-h = || • ||-i + |1 • ||t- 

Finally, let us now extend the operator L defined on C by (12). For any A > 0, consider the continuous 
bilinear form on H x H that coincides on C x C with 

Vv,t/> G C, Bx{ip,^P) = A(v,t/>)2 + [Ta + Th]((-^)1/V, (-5)1/^)- 
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Thanks to Assumption 2.3 and the antisymmetry of H, this form is clearly coercive. Thus it defines a 
strongly continuous resolvent operator and consequently, a unique generator L associated to this resolvent 
operator. More precisely, cp £M belongs to Dom(L) if and only if B\{(p,-) is L^-continuous. In this case, 
there exists f e L^{fl) such that B\{cp, •) ~ (f, •)2 and L(p is equal to f — X(p. It can be proved that this 
definition is independent of A > (see [11], Chapter 1, Section 2, for further details). Let us additionally 
mention that the adjoint operator L* of L in L^{Q, it) can also be described through B\. Indeed, Dom(L*) = 
{ip e W,B\{-,ip) is L^(J7)-continuous}. If £ Dom(L*), there exists f £ L'^{n) such that B\{-,ip) = (f, •)2 
and L*</5 is equal to f — Xip. 



Remark 4-1- For each function <p £ C C H, the application litp can be viewed as a function o/IHI_i. Indeed, 
y-ipeC, (Lv5,'f/')2 = -[Ta + T//]((-S)i/V,(-S)^/V) < [^/ + Cf]||¥3||i||^||i. Hence, the application ip 
Litp £ IElI_i can be extended to the whole space M so that, for each function u £ H, we can define Lu as an 
element of M._i even if Dom(L). The same properties hold for y> £ H i-^ Sip £ IHI_i and <^ £ H i-^ Sip £ 
H_i. 



4-2. Existence of a solution 



Let us now investigate the solvability of (19). Let us consider fixed parameters 6,6 >0 and A > and 
introduce the bilinear form ^ on F x F that coincides on C x C with 

Bls{ip,,p) = X{cp,i,)2 + (^0 ([a + H]2?,v,2?,tA)2 - 0{DtV,^)2 + AV')2. (27) 

In what follows, the parameter S is omitted each time that it is equal to 0. So the form q is simply 
denoted by B^. 

The main result of this section is Proposition 4.2, whose proof can be found in [18], Proposition 5.4, 
without modification (except the value of 9, which does not play a part). Let us however sum up the 
strategy. As explained in [18], the difficulty lies in the time degeneracy of the Dirichlet form i?^, which 
satisfies no sector condition (even weak) . To face this degeneracy, we add a viscosity parameter 6 > and 
consider the form B^g, which satisfies a weak sector condition. Thus, with the help of the Lax-Milgram 
theorem, it defines a resolvent operator that permits to solve (ua,5 is the unknown) 

S 



Xux,5 - (L + eDt)ux,5 - \^-jD^ux,s = h. 

It then remains to let the parameter S tend to 0. For this purpose, we have to establish estimates for the 
family ([-DtUA. 5(2)5 as S goes to 0. This can be done when the right-hand side h is regular enough with 
respect to the time variable by using Assumption 2.3. Thus we state 

Proposition 4.2. Suppose that h £ L'^{f2) n Dom(£'t) and d eH. Then, for any fixed 9 >0 and A > 0, 
there exists a unique solution Ua £ F of the equation Xu\ — Lua — 9Dt\i\ = h + d, in the sense that \/ip £ F, 
_B^(ua,¥') — (h,<^)2 + (d,t^)_i^i. Moreover, DtU\ £ H and we are provided with the following estimates, 
which only involve the constants m,C2,C2 of Assumption 2.3 (in particular, they depend neither on X nor 
on 9), 

A|uA|^ + m||uAi|?<M + M^, (28a) 
A m 

A| AuaI^ + mil AuAil? < ^ + 2Mi + 2(C," + C^ f f ^ + ^) / m\ (28b) 

A TO \ X m I I 



In the case d £ L'^{fi), then u\ £ Dom(L). 
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Finally, Ux is the strong limit in H as 6 goes to of the sequence {u\,s)x.5, where U\,s is the unique 
solution of the equation: Vip E F, g{u\ s, (p) = (h, ip)2 + (d, i, and the family {DtUx.s)s is hounded 
inL^{n). 

4-.3. Control of the solution 

We now aim at determining the asymptotic behavior, as A goes to 0, of the solution (in the sense of 
Proposition 4.2) of the equation 

Xu'x-{L + 9Dt)u\ = h,. (29) 

Since A is no more fixed, 6 = 0(A) is not fixed either. More precisely, we are interested in two specific behaviors 
of the function 9 = 0{X). We will focus on both possibilities limA^o d{X) = and limA^o ^(A) — +oo, which 
respectively correspond to a small and big time/space evolution ratio. In both cases, we will show that 
A|u5^|2 and that (V'^u^)a converges in {L'^{fi)Y as A goes to 0. 

Proposition 4.3. Let ihx)x>o be a family of functions in M^iC] L^{f2) which is strongly convergent in M^i 
to bo and which is bounded in Ti. (see definition (26)). Suppose either lim^^o ^(A) = or limA^o ^(A) = +oo. 
Then the solution Ua G F o/ the equation Xxix — Lua — ODtUx = (in the sense of Proposition 4-2) satisfies: 

• there exists r; € D such that (— S)"'^/'^ua rj as X goes to in D, 

• A|ua|2 ^ as X goes to 0. 

The limit r; G D does not depend on the function 9 but only on its limit as X goes to 0. 

The first step of the proof consists in investigating the case when the operator L is replaced by S in (29). 
This situation is more convenient because of the common spectral decomposition of S and Dt . 

Proposition 4.4. Suppose either liniA^o'^(A) = or limA^o^('^) = +oo. Let (bA)A>o be a family of func- 
tions in EI_i that is strongly convergent to bo in H_i. Let (va)a>o be a family of functions in F that solves 
the equation (for any X> 0) Ava — Sva — 9Dt'Vx = bA in the following sense, 

VcpGF, A(vA,C^)2 + (vA,¥')l-e(AvA,'/5)2 = (bA,C^)2. (30) 

Then there exists G D such that A|va|2 ^ and |(— S)^/^va — T7I2 as X goes to 0. The limit 77 G B does 
not depend on the function 9 but only on its limit as X goes to 0. 

Proof. In what follows, the parameter X of 9{X) is sometimes omitted when there is no possible confusion, 
but keep in mind that 9 does depend on X. From Lemmas 5.10 and 5.11 in [18], we can assume that, for 
any A > 0, bA G L'^{Q) n Dom(L)t) n H_i and converges to bo G H-i. From Proposition 4.2. we can suppose 
that va G Dom(S). Remind that — S = /jj^^^x'i?(da;,d?;) and —Dt — j^^^^iyE{dx,dy). Choosing <^ = Va 
in (30); we have 

a|va|2 + !|va||? = (bA, va)2 < qivAlli < c\ (31) 

where C = supA>o l|t)A||-i- Thus we can find 77 G B and a subsequence, still denoted by (va)a, such that 
((— S)^/^va)a converges weakly in L^(J7) to r]. 

Now we claim supA>o ||Ava||-i < 00 and sup;^>o ||^^tVA||-i < 00. 



|(Ava,¥3)2| 



X{X + x + i9y) ^ dE^^^^{Ax,dy) 



<(/ ^^ ^[(A + x)2 + g2^2] di;b.,b.(dx,d,)) (/ xAE^.^{dx,Ay) 



684 R. Rhodes 

1/2 



A>0\JR+xR / 



A>0 \Jr+ X 

= sup||bA||-il|v?||i. 

A>0 

Since 9DtV\ = Ava — Sva — and ||Sva||-i < ||va||i (cf. Remark 4.1), then DtVx e IHI_i and 
supa>o ll^^-DtVAll -1 < oo. So there exists a bounded family (Fa)a>oo of continuous Hnear forms on E) C L^(i7) 
such that VA > 0, V95 G C,Fa((-S)1/V) = e(A)(AvA, ¥')2- Moreover, from (31), (Ava)a converges to in 
L'^{n) so that, yipeC 

Fa((-S)1/V) = (Ava,c^)2 + ((-S)1/2va, (-S)1/V)2 - (bA,<^)-i.i 

(32) 

^(77,(-S)i/V)2-(bo,<p)-i,i 

as A goes to 0. So (Fa)a>o is weakly convergent in D* (topological dual of D) to a limit denoted by Fg as 
A goes to 0. 

We now aim at proving Fo(t7) = 0. For A, /i > 0, using the antisymmetry of the operator Dt we obtain 
Fa((-S)V2v^) = 0(A)(Ava,v^)2 = -0(A)(Av^, va)2 = -^F^((-S)V2v,). (33) 

If limA^o ^(A) = 0, we pass to the limit as A goes to and we deduce Fo((— S)^/^v^) = 0. It just remains to 
pass to the limit as /x goes to and it comes out Fo{t]) = 0. Otherwise, if hmA— o^(A) = +00, we first pass 
to the limit as /U goes to 0, then as A goes to 0, and we also obtain Fo(r;) = 0. 

Let us now establish the limit equation, which connects Fq, r] and bo. From (32), for any ip dC 

(r,, (-S)i/ V)2 - Fo((-S)i/ V) = (bo, (34) 

This limit equation permits us to investigate uniqueness of the weak limit and, as a bypass, the last state- 
ment of Proposition 4.4. Indeed, let 9 and 9' be two functions having the same limit, say limA^o ^(A) = 
limA-,0 ^''(A) = (the case limA-to ^(A) = liniA^o ^'(A) = +00 is quite similar). Consider the two families of 
associated solutions (va)a and (v^)a of Eq. (30), and two possible weak limits h and rj' of respectively two 
subsequences of (va)a and (v^)a. Define the corresponding linear forms (Fa)a and (F';^)a as well as their 
limits Fq and Fg as described above. Then (34) provides us with the following equality: 

V(-S)i/VeD, (r,-r,',(-S)i/V)-[Fo-F[,]((-S)i/V). (35) 
As we proceeded for (33), we can establish 

Fa((-S)V2v;) = -^|^f;((-S)V2va). 

Now we pass to the limit as A goes to along the first subsequence and then as /i goes to along the second 
subsequence. We obtain Fg(j7') = 0. Reversing the roles of A and fi, we also establish Fg(r7) = 0. Choose now 
(-S)i/V = r]-T]' in (35), this yields 

|r,-77'l2 = [Fo-F[,](r,-,7') = 0. 

In particular, the weak convergence holds for the whole family ((— S)^/^va)a towards rj. Let us now tackle 
the strong convergence of (va)a- Choosing <^ = va in (34), passing to the limit a A goes to and using 
Fg(h) = 0, this yields 

(h,h)2 = lim(bo,VA)-i,i = lim(bA,VA)-i,i ^ = ^ lim[A|vA|^ + |1va||?]. (36) 

A — A — *0 A — *0 
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In particular, |h|2 = liniA^o |(~S)^/^va|2- Thus, the convergence of the norms implies the strong convergence 
of the sequence ((— S)^/^va)a to rj in L^(i7). As a bypass, (36) also implies the convergence of (A|va|2)a to 
0. □ 



Outline of the proof of Proposition 4.3. Let us now briefly explain how to deduce Proposition 4.3 from 
Proposition 4.4. At first, note that = (A — L — Dt)~^ (hx). The plan of attack then consists in defining 
the operator Px-.H^H hy Px{h) = (L - S)(A - S - A)"Hb) and in noting that 

(A - L - Dt)-^ = (A - S - Dt)-^[l - Px]-\ 

If II^aIIw^w < 1, then the operator [I — Pa] is invertible and [I— PA]^^(bA) converges inH_i. Thus Propo- 
sition 4.4 applies. Actually, ||Pa||w— ^ 1 but this norm is finite. An iteration procedure shows that the 
general situation boils down to the particular case ||Pa||-h^w < 1- Rigorous details can be found in [18], 
Proposition 5.12. □ 



5. Ito's formula 



Since the generator of the process is not regular enough to work on, we introduce regular approximations 
through viscosity methods. This allows us to get rid of the time degeneracy. Let us consider a standard 1- 
dimensional Brownian motion {B'^-.t > 0} that is independent of {Bt]t > 0} in such a way that {{B'^, Bt);t > 
0} is a standard (d+ l)-dimensional Brownian motion. Let us then define the (d + l)-dimensional diffusion 
process X"^'^, starting from 0, as the solution of the SDE (we still use the convention that {x,X) stand for 
a. {d+ l)-dimensional vector where .t e M and X e M'') 

xf = j\e'^-",b{X:\^))dr+ (^0,j\{xf,^)dBr^+iVSB'„Q). (37) 

The associated diffusion in random medium V^'^ defined by Y^'^{u}) — t—s.soj is a 12- valued Markov process, 

^ t 

which admits tt as invariant measure (similar to Section 3). So it defines a contraction semigroup on L^{f2, tt). 
The associated (non-symmetric) Dirichlet form is given by (27) (with = e^/^"") with domain F x F and 
satisfies a weak sector condition (see [11], Chapter 1, Section 2, for the definition). The generator L'^''^ is 
defined on Dom(L'^''') = {u e F; B^ ^{u, •) is L^(]7)-continuous} (see [11], Chapter 1, Section 2, for further 
details). It coincides on C with L -I- 9Dt + {d/2)D^. Since b and a are globally Lipschitz (Assumption 2.2), 
classical tools of SDE theory ensure that 

Eo 



i2 



sup \{s'^-%xl)-Xt-'f 



0<t<T 



dTT as (5 goes to 0, (38) 



where both diffusions start from 0. The additional time regularity of this setting enables us to apply the Ito 
formula for suitable functions in the domain of L''*, typically functions ux^s given by Proposition 4.2. It 
remains to make the parameter 6 vanish. Technical details are explained in [18]. 



Theorem 5.1. Let f e L'^{fi) and a family (ua)a>o in ^ such that: 

(1) V<peF, i?»(uA,¥') = (f,<p)2, 

(2) for each A > 0, there exists a sequence {vlx,5)&>o ^ that converges in M towards Ux - Moreover (uA,d-)5>o 
satisfies Bl g{ux.s, •) = (f, ■)2, 

(3) for each fixed A > 0, (DtUx,s)s is hounded in L^(f2), 

(4) each function Ux.s has continuous trajectories, that is, for /i almost every ui ^ Q , the function (t, x) G 

^ \ix.s{Tt,x'^) is continuous. 
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Then, choosing 9 (X) ~ "/(2/3) and X = e'^^ , P-^ ci.s., the following formula holds 



f)(y;)dr + 



v"u*.,(y/)dB, 







where Ptt is i/ie /aw of the process starting with initial distribution n on f2. 
6. Ergodic theorem 

Wc now focus on the ergodic properties of the family of processes {Y^)^. Basically, the difficulty lies here in 
the difference of evolution rates of the time and space variables. The strategy will consist in establishing an 
orthogonal decomposition of LF'{{T) that allows, in a way, to separate the variables and to exploit separately 
their ergodic properties. Concerning the space variables, this is carried through from Assumption 2.4. 

Theorem 6.1. LetieL^{n). Then 



Proof. Decomposition of the space L'^{Q). Let us first establish the following orthogonal decomposition of 
the space L/^^fl) 

L^{Q) = Inv e Closure(R), 

where R {Sy?; ip G Dom(S)} n Dom(A) and Inv = {f S L'^{Q)-yx e R'*, To.^^f = f}. It is readily seen that 
Inv and R are orthogonal and that L'^{Q) D Inv ® Closure(R). 

Conversely, let f £ L'^{Q) such that VSv? G R, (f , ^ip)2 = 0. Consider h G L'^{f2) n Dom(Di). Hence Propo- 
sition 4.2 provides us with a family (ua)a G F n Dom(S) such that B^{u\, ip) ~ (h, ip)2 for any G F. We 
emphasize that the equality 

Aua - Sua = h (39) 

also holds in the sense, so that Sua = Aua — h G R. From (28a), (Aua)a is bounded in L'^{fi). Consider 
a weak limit h* G L'^{Q). Let us prove that Sh* = 0. Let ip G Dom(S). Integrate (39) against ip, multiply 
by A, pass to the limit as A goes to and obtain (h*,S<p)2 =0. Since S is self-adjoint, h* G Dom(S) and 
Sh* = 0. From Assumption 2.4, h* G Inv. Let us additionally mention that h* is the orthogonal projection 
of h into Inv. Indeed, integrate (39) against G Inv, note that (Sua, v)2 — 0, pass to the limit as A 
goes to and obtain (h — h*,(p)2 = 0. In particular, (Tg.^h)* ~ h* for x G K''. Returning to our first 
objective, observe that (f, h)2 = (f,AuA — Sua)2 = (f,AuA)2. Passing to the limit as A goes to 0, we obtain 
(f, h)2 = (f , h*)2 for each function h G L2(J7) nDom(£>(). We are now in position to conclude: for each a; G M'^ 
and h G L'^{Q) n Dom(Di), 



Since L^{il) D Dom(£'t) is dense in L^{f2), we deduce To^xi = f • The decomposition follows. 

Particular case. Let us first prove Theorem 6.1 in the case f = f * + S<p, where <p G Dom(S), S<p G Dom(Dt) 
and f* G Inv. For each A > 0, Proposition 4.2 provides us with a family (fA)A that satisfy (in the sense of 
Proposition 4.2) 




f (y/) dr - tTr{f) ^0 as e goes to 0. 



(To.,f,h)2 = (f,To,_,h)2 = (f,(ro,_,h)*)2 = (f,h*)2 = (f,h)2. 



AfA-LfA-0(A)AfA = Sc^. 



(40) 
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Choose now A = e^^ and e{X) = X^-^/C^P), Theorem 5.1 yields 



t/e ' t/e ' 

Jo Jo 



in such a way that, using the decomposition f = f * + Sip, 



,2/3 



t/e" 



tie- 



2/3 



f(y;)dr = e"y r(T^,ow)dr + e2/5 

Jo 



r-tle^l' 



Let us now prove that the above right-hand side behaves as its first term as e goes to 0. From Proposition 
4.3 (note that Sy G IHI_i), we have A^|fA|2 + A||fA|ii — > as A goes to 0. Thanks to the Jensen inequahty, the 
invariance of the measure tt and Proposition 4.3 (note that Sy? G H-i), we deduce 



-2/3 



t/e ■ 



Similarly, E^[(e2^(f^2^f(Fj^^20) — fe^^f (tj)))^] < 2e^^\^^'2ii\^ and this latter quantity converges to as e goes to 
0. Concerning the martingale part, we just have to estimate its quadratic variations 



,2/3 



/ V"f,24K;)dB, 
Jo 













e^^ f 






Jo 



/ |v-f,2^(r;)|'dr 

Jo 



to conclude that it also converges to 0. Thus it remains to study the convergence of the first term. The 
classical ergodic theory ensures that £"/o^^ f*(7r,ow)dr converges in L'^{f2,TT) to a limit function that is 
invariant under time translations. Since f* is already invariant under space translations, it is readily seen 
that this limit function is also invariant under space translations and consequently constant (see Definition 
2.1). Thus it is equal to 7r(r) = 7r(f). 

General case. Consider now a function f G L^{f2) and denote by | • |i the i^(i7)-norm. Obviously, we can 
find a sequence (f„)„ C L^{Q) that converges in L^{f2) towards f. Then, using the invariance of the measure 

TT, 



E, 



,2/3 



<E, 



pt/e - 

/ f(y/)dr-t7r(f) 
Jo 

/ [f-f„](r/)dr 
Jo 

ft/e^" 

■^n f„(r/)dr-t^(f„) 
Jo 



■ t|7r(f) - 7r(f„ 



<2i|f-f„|i+E, 



-2/3 



t/e' 



f„(r;)dr-i7r(f„) 



It just remains to choose n large enough to make the former term in the above right-hand side small and 
then to choose e small enough to treat the latter term. This completes the proof. □ 
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7. Invariance principle 

Notations. From now on, fix 6{X) = X^^'^/'^'^P) . For i e {1, . . . , d} and A > 0, Proposition 4.2 provides us with 
a solution e F n Doin(L) of the equation 

Au5,-L<-0(A)Aul=b„ 

and with a sohition ca £ F n Dom(L) of the equation 

Aca - Lca - 0{X)DtCx = c. 

An integration by parts proves that b,c fulfill the assumptions of Proposition 4.3 (see [18], Lemma 5.14, 
for the detailed proof). Thus it makes sense to define = limA^o V'^u\ and k = limA^o V'^ca, where both 
limits are understood in the L'^{Q) sense. Moreover, from Proposition 4.3, A|ua|2 + A|ca|2 as A tends to 
0. 

Consider now a bounded continuous function / e C{W^) and the solution z^^^ of Eq. (2). The Feynman- 
Kac formula gives a probabilistic representation of z^ ^j (cf. [17], Theorem 3.2 and Remark 3.3) 

z,^^{t,x)^¥.,[f{XI)cMQl)l (41) 

where Qf = j^[e^^c + d]{r /e°' ,X^/e^ ,uj)dr. As guessed by the reader, the strategy consists in studying 
the convergence in law of the couple of processes {X^,Q^) and then in establishing a uniform integrability 
argument. 

Let us first tackle the convergence in law of the couple {X\Q^). Define Q^ = y^[c + ePd]{Y^)Ar and 

remind, from Section 3, that (X^,Q'^) and {e^ X^,^^2i3 ,e^Q^.fg2fi) have the same law, where both processes X^ 

and X^ start from 0. 

Applying Theorem 5.1 to the functions u^2fi and c^2f3 yields 

ePXt/^2, = HI + / (cr + V^u*., )(y/) dB„ 
Jo Jo 

where 

ilf =£3/3 / u,2,(y/)dr-eV^4>;%2^)+e''u,2,H, 
Jo 

Gl = e^" / c,2, (r/) dr - 61^0,2, ) + e'''c,2, (w). 

Jo 

For the sake of clarity, it is worth recalling that Pjr denotes the law of the process starting with tt as 
initial distribution and that the measure tt is invariant for the process . Let us now establish that the 
finite dimensional distributions of both processes H'^ and converge in P7r-probability to 0. Using the 
Cauchy-Schwarz inequality and the invariance of the measure tt, we obtain 

E^[iHlf]<3{2 + t')s^P\u,2,\l 

and this latter quantity converges to as e goes to 0. Likewise, E7r[(G'f )^] converges to as e goes to 0. Then, 
from Theorem 6.1, the process e'^^ d(y/) dr converges, at least in probability, to the deterministic pro- 

cess 1 1~-> t7r(d). Finally, we tackle the convergence of the martingale part of the process {£^X"./^2i3,e^Q'^./^2i3), 
which matches 

/ [(a + V'^u:.,),V%%,](y/)dE,. 
Jo 
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Clearly, under P^, the difference between this latter process and £^ Jq + '^*](^/) di3r vanishes as 

e goes to in P^r-quadratic mean. The c^uadratic variations are easily computed 



,2/3 



t/e 



(y/)dr. 



Using Theorem 6.1 again, these quadratic variations converge, at least in Pjr -probability, towards the deter- 
ministic process 1 1-^ At, where the nonnegative symmetric matrix A is given by 



A: 



K*(cr + ^*)*(w) k*k{uj) 



d7r(ti;) 



To sum up, the finite dimensional distributions of the process {e^ Xl^^2i3 , converge in law to the 
1/2 

process (0, D)t + A B'^., where D = 7r(d) and B' is a (d -I- l)-dimensional Brownian motion. Actually, this 
convergence holds in the sense of weak convergence of processes in the space C([0, r];R''+^) for each fixed 
T > 0. Section 8 is devoted to the proof of this fact. 



8. Tightness 



We keep the notations of Section 5. To prove the tightness of the underlying stochastic processes, we use 
Garsia, Rodemich and Rumsey inequality (GRR's inequality) (cf. [19]). More precisely, we follow [14], Section 
3. Nevertheless, the required condition (43) is stronger than the one stated in [18], Section 9, that is g S H_i. 
But this methods provides us with a uniform integrability criterion needed in Section 9. Once again, we 
work on the viscosity approximations Y^'^ to avoid facing the time degeneracy of the generator of Y'^ . 
From [14], Theorem 3.2, for a given function g S L°°{n), the function 



rig)it,Lo) = E 



exp(^^ g(F;'')dr 



belongs to Dom(L'^*) and satisfies dtr{g) = L^'''_r(g) + g-r(g) with initial condition r{g){Q,uj) = 1. Then 
[14], Proposition 3.3, ensures that 



|r(g)(t,.)l2<cxp(2tSp(L^''* + g)), 



(42) 



where Sp(L^''' + g) = sup|^|2^i(v, [L'^-'^ -t- g]¥')2 and the sup is taken over (p G Dom(L'^''') C F. Thus, for any 

, we have 

{Yr''') dr) = riae^g){e-^P{t - s), •) 



a > 0, by the invariance of the measure tt, we have 



exp I ae 



r/3 



<|r(ae'^g)(e-^'^(t-s),.)l2 
< exp((t - s)Sp(e-2''L-'^ + 



ae 



It just remains to give a bound for Sp(e ^^L'^'"' + ae '^^g) and we want it to be independent of 5 and e. 
Suppose now that 



V(peC, |¥:'|2 = 1 =^ (g,¥'^)2<C(43)||¥'lli 
for some constant (^^43) > 0. For any e F with \ip\2 = 1, 



(43) 
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< 



-2/3 



m\\ip\\l + ae '^^g^Hv 



< 



4m 



Let us additionally assume that, for each fixed lo, the mapping {t,x) i— > g{t,x,u)) is globally Lipschitz. From 
(38), we can then pass to the limit as 5 goes to in the inequality 



exp ae 



r/3 



and get the bound 



exp I ae 



r/3 



< exp 



(t-s)a2c2 



(43) 



Am 



;(r/)dr 



< 



exp 



{1-3)0^0^ 



(43) 



4to 



(44) 



Let us now state GRR's inequality, whose proof can be found in [19] or [14], Proposition 3.1: 



Proposition 8.1 (Garsia— Rodemich— Rumsey's inequality). Let p and W he strictly increasing con- 
tinuous functions on [0,+oo[ satisfying p{0) = !f'(0) = and limt^oo ^'(^) = +00. For given T > and 
f G C([0,r];K''), suppose that there exists a finite B such that; 



T pT 



Asdt < B <oo. 



m-g{s)\ 
lo Jo ^ V Pi\t-s\) 

Then, for all < s <t <T, 



\g{t)-g{s)\<8 I ^^-^f^)dp{u). 



(45) 



(46) 



Choose now W{t) = e* - 1, ^~^{t) = ln(l + 1) and p{t) = ^/t. As explained in [14], Section 3, (46) provides 
us with the following estimate of the continuity modulus: 



sup \g{t)-g{s)\<8VS\niS-') 

\t~s\<S 
0<s<t<T 

r/3 r*/= 



In 4 



T i-T 



^0 



exp 



W)-9{^)\ 



dsd< +6 



(47) 



Choosing g{t) = g(^/) dr, taking the expectation and using Jensen's inequality, we obtain 











sup 


en g(r/)dr 




|t-s|<5 






0<s<t<T 





In 4 



T i-T 



E. 



exp 



P/,XCg(n^)dr| 



dsdt +6 



Gathering (44) with the above estimate, we finally obtain the desired continuity modulus estimate 



E, 


sup 






|t-s|<5 





pt/e - 

/ g(n')dr 



(48) 



0<s<t<T 
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for some positive constant Ct- Note that both functions b and c satisfy (43). Indeed, for any cp E C (the 
case g = b is similar) 

= -2j2ivii,^vD,v)2<MmL-in)\v\2Mi. 

In particular, (48) holds for g = b and g = c. Note that, from the boundedness of d and cr, the 
Burkholder-Davis-Gundy and the Kolmogorov criterion, the tightness of the processes e^'^/g d(y/)dr 

/ 2/3 

and ef^J^'" cr(y/)dB^ raises no particular difhculty. As a consequence, the tightness of the processes 
e^X^,l^2i3 and e^C^.j^ifi is proved. 



9. Girsanov's transform and limit equation 



To determine the limit of ^^ [^(i, x) (see (41)), we now aim at applying the convergence in law of the 
couple (e^X^i^iii , e^Cfij^nii ) to the function g{x\^ . . . , Xd+i) = f{xi, . . . , Xd)e^''^^ for some continuous bounded 
function / € Cb(]R'';M). This requires a uniform integrability argument, which is derived from (44), applied 
with g = c. Since d is bounded, the random variable exp(£^'^ /^^^ d(y/) dr) is bounded too. So the random 
variable /(e'^Xj /^2/3) exp(e'^Q( 7^2/3) is uniformly integrable and thus converges in 7r-probability (cf. Section 



1/2 

7) to Eo[g((0,D)t + A B[)]. We now aim at finding a PDE that characterizes this limit. Let us consider 



the 



-de-composition in blocs of A and A 



1/2 



A 



An A 



12 



A21 I A 



22 



(A''\2 



By the Girsanov transform, we can define a new probability P on C {[Q ,T];W^'^^) as follows 



= exp (A^f -tJ,')B[~M2\ 



— — \I2 1/2 

where Tt denotes the natural filtration of the Brownian motion. Under P, Bt — B'^. — {A^2 1^22 )*^ i 

1/2 ~ 

Brownian motion. Let us now rewrite Eo[g((0, D)t -\- A B'^)] in terms of B 

M9{{Q.D)t + A"^B[)] 

= ni{(A\'i\A]'2)B't)cMDt + (^2^,^2^)50] 



IS a 



where 



. , ,-^1/2 -rl/2, ~ ,-rl/2-;-l/2 -rl/2-^1/2,, —rt 

/((Vi .Ai2 )Bt + t{A^i A2 +^12 ^22 ))cx^l>{Dt + A22- 



= E[/((4f 4f )i?,+Ct)]e^*, 



C= / {(T + C)i^<^T^ and U = 



d dTT. 



(49) 
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We point out that the quadratic variations of {Aj^^ Ay2, )-Bt are equal to 

A = An^ f (^ + r)(T + r)*d7r. (50) 

Hence there exists a standard d-dimcnsional Brownian motion Bt such that 

Eo[gmD)t + A'^^B[)]^Eo[fiCt + A'/^Bt)c^% 
To sum up, we identified the limit in tt probability of 

z,,UQ,t)=Eo[f{X!)exp{QI)]=Eo[f{e^xl/,.,)exp{e^Ql/,.,)] 

asEo[/(Ct + Ai/2B,)c^*]. 

Let us now to determine this limit when the starting point is not but x e M'': 

E4f{XnoMQn] ^'=™?^^3.2 Eo[/(x + Xr'<"-^^''>")exp(Qr<«-^'^>'^)] 

in law with^cspect to ^ ^^^^^^ ^ exp(Qf''^)] 

TTprob Eo[f{x + Ct + A^/'^Bt)c^% 

To complete the proof of Theorem 2.6, it just remains to explain why the coefficients A, C, U only depend 
on the 3 cases a — 2/3 < 0, a — 2/3 > or a — 2/3 = 0. Since ^ and k only depend on the case a — 2/3 < 0, 
a — 2/3 > and a — 2(3 — (see Proposition 4.3 for the first two cases and [18] for the last one), the same 
property holds for A, C and U, which can be expressed in terms of ^ and k (see (49) and (50)). As a 
consequence, this completes the proof of Theorem 2.6. 
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